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Abstract 

The detectability of a possible “high” energy (> 100 GeV) component of 
Gamma-Ray Bursts (GRBs) using AMANDA/ICECUBE large area muon de¬ 
tectors is examined within the context of a specific model of such high energy 
gamma ray production within GRBs, namely, the proton-synchrotron model, 
which requires protons to be accelerated to ultrahigh energies > 10 20 eV within 
GRBs. In this model, the high energy component is distinct from, but may well 
be emitted in coincidence with, the usual “low” (keV-MeV) energy component 
observed by satellite-borne detectors. The AMANDA/ICECUBE detectors can 
detect TeV photons by detecting the secondary muons created by the TeV pho¬ 
tons in the Earth’s atmosphere. We calculate the muon signal to noise ratio 
in these detectors due to TeV gamma-rays from individual GRBs for various 
assumptions on their luminosity, distance (redshift), Lorentz Gamma factor of 
the underlying fireball model, and various spectral characteristics of the GRBs, 
including the effect of the absorption of TeV photons within the GRB as well as 
in the intergalactic infrared radiation background. The intergalactic absorption 
effect essentially precludes detection of TeV photons in the AMANDA detector 
for reasonable values of the luminosity in the high energy component, but they 
may well be detectable in the proposed ICECUBE detector which may have 
an effective area for downward-going muons a factor of 100 larger than that 
in AMANDA. However, even in ICECUBE, only relatively close-by GRBs at 
redshifts < 0.05 or so can be expected to be detectable with any reasonable 
degree of confidence. We discuss the requirement on the luminosity of the GRB 
in the high energy component for its detectability in ICECUBE. 
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1 Introduction 


Gamma Ray Bursts (GRBs), short and intense bursts of photons observed mostly in 
the energy range of few tens of keV to few MeV with burst duration lasting from 
fraction of a second to several hundreds of seconds (see e.g., Ref. [1] for a review), 
are one of the most powerful and enigmatic astrophysical phenomena in the Universe. 
An afterglow phase, lasting from several hours to several days, has also been detected 
for many bursts. Detections of these afterglow phases in optical and/or x-rays and 
in some cases in radio, have provided valuable information regarding the nature of 
GRBs [2]. For example, large redshifts of the host galaxies of some GRBs, measured 
from observations of their optical afterglows, have unambiguously established that 
GRB sources are extragalactic and indeed lie at cosmological distances. 

The total energy emitted in a GRB, estimated from the observed fluence and mea¬ 
sured redshift (distance) of the burst source assuming isotropic emission, varies over 
a wide range from burst to burst, from about few xl0 50 ergs to 10 54 ergs. On the 
other hand, the material emitting the radiation could itself be highly collimated in 
the form of a jet, making the emitted radiation highly non-isotropic. One possible 
diagnostic, albeit not always a robust one, of a collimated emission is an expected 
break feature [3] in the light curve of the optical afterglow of the GRB. Collimated 
emission significantly reduces the estimated total energy of the burst compared to 
that estimated assuming isotropic emission. From an analysis of the measured red- 
shift and estimated collimation angle (from the break feature in the optical afterglow 
light curve) of some 18 bursts, Ref. [4] finds a rather narrow distribution of the total 
energies of bursts with typical values around 5 x 10 50 ergs with a FWHM of a factor 
of ~ 5. 

The satellite-borne detectors, such as the BATSE on board the Compton Gamma Ray 
Observatory (CGRO) and earlier detectors on board the Vela series of satellites, have 
detected GRBs mostly in the sub-MeV energy region. However, several GRBs with 
emission beyond 100 MeV, and in coincidence with the sub-MeV burst, have been 
detected [5, 6, 7] by the EGRET instrument on board the same CGRO, including 
the long-duration burst GRB940217 [6] with emission extending to ~ 18 GeV. This 
suggests the possibility that GRBs may, in addition to the sub-MeV photons, also 
emit much higher energy photons, perhaps even extending to TeV energies as in some 
highly energetic Active Galactic Nuclei. For power-law spectra falling with energy, the 
photon number flux at TeV energies may be too low for these more energetic photons 
from GRBs to be detected by the satellite-borne detectors which have limited sizes. 
However, ground-based detectors can in principle detect TeV photons from GRBs by 
detecting the secondary particles comprising the “air showers” generated by these 
photons in the Earth’s atmosphere. 
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Indeed, three major ground-based gamma ray detectors, the Tibet air shower array [8], 
the HEGRA-AIROBICC Cherenkov array [9] and the Milagro water-Cherenkov de¬ 
tector [10] have independently claimed evidence, albeit not with strong statistical 
significance, for possible TeV y-ray emission from sources in directional and temporal 
coincidence with some GRBs detected by BATSE. The estimated energy in TeV pho¬ 
tons in each case has been found to be about 1 to 2 orders of magnitude larger than the 
corresponding sub-MeV energies measured by BATSE. More recently, a search [11] 
for sub-TeV (> 10 GeV) gamma rays from GRBs with the Project GRAND array 
of muon detectors reported one possible detection at the 2.7cr level in coincidence 
with one BATSE GRB (971110), and an analysis [12] of the energetics of this event 
within the context of various possible models of TeV gamma ray emission from GRBs 
indicates a significantly higher estimated energy in TeV photons (by four orders of 
magnitude or more) than in sub-MeV photons. 

Note that, since TeV photons are efficiently absorbed in the intergalactic infrared (IR) 
background due to pair production [13], only relatively close by (i.e., low redshift) 
GRBs, for which the absorption due to IR background is insignificant, can be observed 
at TeV energies, which would explain the fact that only a few of the BATSE-detected 
GRBs in the fields of view of the individual ground detectors have been claimed to 
be detected at TeV energies. 

The significantly higher total energies estimated for the TeV photon emitting GRBs 
claimed to be detected by the ground-based detectors mentioned above raise the 
possibility that the TeV photons may constitute a separate ‘high’ (GeV-TeV) energy 
component of GRB emission distinct from the ‘low’ (keV-MeV) energy component 
detected by BATSE. For, otherwise, if the ‘observed’ TeV photons came from the 
same low energy spectrum continued to TeV energies, then it would be difficult to 
understand the significantly higher total energy estimated for the bursts from the TeV 
observations than that estimated for the same bursts from sub-MeV observations by 
BATSE detector, unless the spectrum is hard with a differential power-law spectrum, 
dN/dE oc E~ a , with a < 2, which is not the case for the observed bursts. 

Confirmation of TeV gamma ray emission from GRBs would have major implica¬ 
tions not only for the physics and astrophysics of GRBs as such, but also for several 
other phenomena. For example, some of the proposed mechanisms [14, 15, 16, 17] 
of possible TeV photon emission from GRBs discussed below are directly linked with 
the proposed scenario of GRBs being possible sources of Ultrahigh Energy Cosmic 
Rays (UHECR) [18]. It has also been suggested [19, 20] that while the TeV photons 
emitted by GRBs at large redshifts would be absorbed through e + e _ pair production 
on the intergalactic infrared background, the resulting electromagnetic cascades ini¬ 
tiated by the produced pairs could produce the observed extragalactic diffuse gamma 
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ray background in the GeV energy region. It is, therefore, important to study the 
detectability of the possible TeV photon emissions from GRBs in as many different 
ways as possible. 

In this paper, we study the detectability of possible TeV photons from GRBs in the 
proposed ICECUBE [21] class underground muon detectors, within the context of 
a specific model of such high energy gamma ray production within GRBs, namely, 
the proton-synchrotron model discussed in detail below. The ICECUBE type muon 
detectors can detect TeV photons through detection of the secondary muons that 
are produced in Earth’s atmosphere by high energy photons of energy above a few 
hundred GeV. 

The ICECUBE is the proposed full-scale version of the already operating 
AMANDA [22] detector. These are facilities primarily for detecting high energy 
neutrinos from cosmic sources. Located under the antarctic ice sheet, and consisting 
of arrays of photomultiplier tubes (PMTs) attached to strings and frozen in the deep 
ice, these detectors detect high energy neutrinos by detecting the muons produced 
by the neutrinos within the detector volume in the ice; the muons themselves are 
detected and their tracks (and thereby the arrival direction of the parent neutrinos) 
reconstructed using the array of PMTs which register the cherenkov light emitted 
by the muons traveling through the ice. In the final analysis, however, these are 
just muon detectors, and can, in principle, also detect the muons produced by TeV 
photons in the Earth’s atmosphere. There is one basic difference, though, in the 
detection method employed for photons as against that for neutrinos: While for neu¬ 
trinos one looks for upward-going muons in the detector, photons have to be detected 
only through detection of the downward-going muons in the detector. 

The advantages of using these muon detectors over the conventional air-Cherenkov 
telescopes for doing TeV gamma ray astronomy have been expounded in Refs. [23, 24], 
Compared to air-Cherenkov telescopes, these detectors typically cover much larger 
fraction of the sky with large duty cycles. For example, the AMANDA detector 
covers more than a quarter of the sky with essentially 100% efficiency. The detector is 
sensitive to muons with energies of a few hundred GeV and thereby to parent photons 
of energy of order a TeV and above. The fact that air-showers created by gamma 
rays are relatively “muon poor” (as compared to air-showers generated by hadrons) 
can be, to a large extent, compensated for by the relatively large effective area of a 
detector such as ICECUBE (~ lkm 2 ). In addition, as pointed out in Ref. [24] (AH, 
hereafter), despite the otherwise large background of down-going cosmic-ray induced 
atmospheric muons, the above muon detectors can in principle detect TeV gamma 
rays from transient sources like GRBs because, with the information on the time and 
duration of a burst provided by satellite observation, the background is significantly 
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reduced because it is integrated only over the relatively short duration of a typical 
GRB which is of order few seconds. 

Of course, possible TeV photons from GRBs can be detected by “conventional” air- 
shower detectors, as exemplified by the claimed detections [ 8 , 9, 10] already mentioned 
above. In this paper we focus on the AMANDA/ICECUBE detectors for the advan¬ 
tages mentioned above and in particular for their potentially unique capability of 
simultaneously being detectors of neutrinos as well as photons of TeV energy and 
above, albeit not from the same source. 

In AH [24], the expected number of muons and the signal-to-noise ratios in AMANDA 
and Lake Baikal detectors due to TeV gamma-rays from individual GRBs were calcu¬ 
lated for various different values of the GRB parameters. AH assumed a single power- 
law energy spectrum for individual GRBs continuing from the MeV BATSE band to 
TeV energies. However, as discussed above, it is more likely that the TeV photons 
constitute a separate high energy component distinct from the low energy “BATSE” 
component. In this paper, therefore, we consider such a model of TeV gamma ray 
emission from GRBs, namely, the proton-synchrotron model, in which, the TeV pho¬ 
tons come from a new component different in origin than the sub-MeV “BATSE” 
photons. We study the detectability of the TeV photons in AMANDA/ICECUBE 
detectors within the context of this model, and study the sensitivity of the signal- 
to-noise ratio to various parameters such as the total luminosity of the GRB in the 
high energy (GeV-TeV) component, the redshift of the burst, etc., and in particular, 
the Lorentz Gamma factor (T) of the ultra-relativistic wind in the underlying fireball 
model (see below) of the GRB. 

The rest of this paper is organized as follows: In Sec. 2, we briefly discuss the moti¬ 
vation for considering the specihc model of TeV photon production in GRBs in this 
paper, namely, the proton synchrotron model. The model is discussed in detail in 
Sec. 3. The internal optical depth of the produced high energy photons due to the 
process 77 —> e + e~ within the GRB environment is calculated in detail in Sec. 4 giv¬ 
ing the spectrum of high energy photons escaping from the GRB source. We also give 
the expression for the spectrum of photons arriving at Earth which includes the effect 
of the intergalactic absorption of TeV photons due to the process 77 —> e + e _ taking 
place on the photons constituting the intergalactic infrared background. In Sec. 5, 
the parametrization of the number of secondary muons in photon-induced showers 
in Earth’s atmosphere used in this paper is discussed, and the muon signal-to-noise 
ratio in a ICECUBE class detector is calculated for individual GRBs as functions of 
various GRB parameters. The main results are discussed in Sec. 6 , and conclusions 
are presented in Sec. 7. 
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2 Mechanisms of TeV Gamma Ray Production 
in GRBs: Motivation for Proton-Synchrotron 
Mechanism 

The fundamental source of the radiation from GRBs is thought to be the dissipation 
of the kinetic energy of ultra-relativistic (Lorentz factor T > few 100) bulk flow of 
matter caused by emission from a central engine, the nature of which is not yet 
known. Ultra-relativistic bulk flow of matter is obtained naturally in the fireball 
model (see. e.g., [25] for reviews), in which the central engine initially releases a large 
amount of thermal energy E th mostly in the form of an electron-positron-photon 
plasma in a small volume in a relatively baryon dilute medium. The baryon load - 
consisting of the baryons emitted by the central engine as well as those present in the 
ambient medium — is assumed to be small, i.e., r)~ l = Me 2 / E t h <C 1, where M is the 
total mass of the baryon load. Because of the high optical depth of the photons due 
to the pair production process 77 —> e + e _ , the photons cannot escape; instead, the 
radiation pressure pushes the baryons into a shell of matter that expands outward. 
The initial thermal energy of the fireball thus gets transferred to the kinetic energy 
of an expanding shell of matter which, because of its small baryonic load, eventually 
attains Lorentz factor T ~ 77 3 > 1 when essentially all of the initial thermal energy of 
the fireball is converted to kinetic energy of bulk flow. If the emission from the central 
engine is unsteady and sporadic, then the above process may result in multiple shells 
of baryons moving with different lorentz factors, whereas a single ultrarelativistic 
baryonic shell is obtained if the emission from the central engine is steady over the 
duration of the emission. In either case, the observed radiation is believed to arise 
from reconversion, through some dissipative process, of the kinetic energy of the ultra- 
relativistic bulk flow of matter back to internal energy at a large radius beyond the 
“Thomson photosphere” where the optical depth due to pair production, Thomson 
scattering and other processes is small. 

The dissipation of kinetic energy of bulk flow may occur, in the case of unsteady 
emission, when different shells of matter moving with different T collide with each 
other and give rise to “internal shocks”, or, in the case of single impulsive emission, 
when the ultrarelativistic shell of matter decelerates as it sweeps up matter from the 
external medium, giving rise to an “external shock”. The multiple baryon shells in 
the case of unsteady emission from the central engine would collide and eventually 
merge into one shell, which pushing into the external medium will also give rise to 
an external shock. Electrons as well as protons can be accelerated to high energies 
through Fermi mechanism at the internal and external shocks. In the multiple shell 
scenario, the prompt GRB phase is thought to result from the synchrotron radiation 
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(and possibly also inverse-Compton scattering) of the electrons accelerated in the 
internal shocks while the afterglow phase results from the same process(es) occurring 
at the external shock. In the single shell scenario, on the other hand, there are no 
internal shocks, and both the prompt and afterglow phases result from processes at 
the external shock as it moves through the external medium. Although currently 
internal-external shock scenario is popular (see Ref. [25] and references therein), it 
has been argued [26] recently that the single impulsive external shock scenario can 
efficiently reproduce the observed properties of the prompt as well the afterglow 
phases of GRBs. 

As mentioned in the previous section, a possible TeV component of GRB radiation is 
most likely to be of a different origin than the process that gives rise to the keV-MeV 
burst component. Mainly three processes have been suggested so far which can give 
rise to TeV photon emission from GRBs within the context of the fireball model of 
GRBs: (1) Synchrotron self-compton (SSC) process (see, e.g., [27]), i.e, the inverse 
compton scattering of sufficiently high energy electrons on the ambient synchrotron 
photons within the source, the synchrotron photons being radiated by the same non- 
thermal electron population; (2) photo-pion production by ultrahigh energy protons 
within the source and the subsequent decay of the neutral pions to photons [28, 29, 
30], and (3) synchrotron radiation by ultrahigh energy protons above 10 20 eV in the 
magnetic held within the source [14, 15, 16, 17]. In general, the required high energy 
non-thermal electrons and protons can be accelerated either in the internal shocks or 
in the external shock, and the associated TeV photons can be emitted either in the 
burst phase or the afterglow phase. In the present paper, we shall consider the case 
of prompt TeV photon emission in coincidence with the burst phase. 

Both photo-pion production and the proton-synchrotron processes are motivated by 
the suggestion [18] that GRBs may be the sources of the observed ultrahigh energy 
cosmic rays (UHECR) [31] whose spectrum extends beyond 10 20 eV. In this scenario, 
UHECR particles are assumed to be protons accelerated within GRBs. Thus, in 
this scenario, while the synchrotron radiation of electrons produces the observed 
low energy (keV-MeV) component, a high energy (GeV-TeV) component could be 
produced by the ultrahigh energy proton component either through synchrotron or 
photo-pion processes. 

If one assumes that the total energy density of the high energy GeV-TeV component 
photons within the source is significantly higher than that of the low energy “BATSE” 
photons, as seems to be the case for the TeV photon emitting GRBs claimed to 
have been detected so far, and if one further assumes that the energy density of the 
magnetic held within the source is roughly in equipartition with the total energy 
density, then it can be shown [12] that, for production of photons of observed energy 
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~ TeV and above, the proton-synchrotron process dominates over the photo-pion 
process, both processes being clue to the same underlying non-thermal proton content 
within the source. On the other hand, if the SSC process were to produce the high 
energy component with a significantly higher energy density content than that of 
the low energy component, both components in this case being due to the same 
underlying non-thermal electron content within the source, then the energy density in 
the magnetic held would have to be significantly less than its equipartition value [12], 
a condition that is at variance with models in which GRBs are sources of UHECR, 
which by and large require a magnetic field within GRB close to the equipartition 
limit [18]. 

In this paper, we shall work within the framework of GRBs as sources of the observed 
UHECR up to ~ 10 20 eV and, based on arguments mentioned above, as far as the 
question of production of a high energy photon component of GRBs is concerned, we 
shall limit ourselves in the rest of this paper specifically to the proton-synchrotron 
model. 

The proton-synchrotron process may also play an important role in the internal en¬ 
ergetics of GRBs. For example, it has been argued (see, e.g., [16, 17, 20, 32]) that, 
if the fundamental source of energy of the GRBs is indeed the kinetic energy of the 
ultrarelativistic bulk flow of matter as in the fireball model, then, at least initially, 
one would expect the total energy content in protons to be higher than that in elec¬ 
trons by a factor of ~ m p /m e ~ 2000, where m p and m e are proton and electron 
rest mass, respectively. Immediately after dissipation of the kinetic energy of the 
bulk flow through formation of internal and external shocks, the total energy con¬ 
tent of protons would, therefore, be higher than that of electrons. Transfer of energy 
from the protons to electrons through coulomb interactions is not efficient [32], In¬ 
stead, as suggested by Totani [32], protons accelerated to energies above 10 20 eV may 
first produce TeV energy synchrotron photons which, in collision with the sub-MeV 
photons within the source, may produce e + e _ pairs through 77 —> e + e~, thereby 
providing another channel through which energy can be transferred from the protons 
to electrons. If this process is very efficient in a GRB, then that GRB will not be a 
significant emitter of TeV photons. If on the other hand, the above pair production 
process is not very efficient, then the TeV photons would escape, resulting in strong 
TeV photon emission. The efficiency of the above pair production process for the 
TeV photons turns out to be a sensitive function of the initial bulk Lorentz factor 
T (see below). Thus, small fluctuations in T from burst to burst could give rise to 
wide variation of the TeV photon output from burst to burst even if the initial total 
energy input in all GRBs lies in a narrow range. Of course, the total energy output 
of a GRB in TeV photons will also depend on the shape of the accelerated protons’ 



spectrum: If the proton spectrum is sufficiently hard so that most of the total energy 
in the proton component lies at the highest energy end of the spectrum where the 
synchrotron emission process is efficient, and if the pair-production process for the 
resulting TeV synchrotron photons is inefficient, then the bulk of the energy in the 
proton component (which is initially higher than that in the electron component) will 
escape from the source in the form of TeV photons, which may naturally explain the 
significantly higher total energy in the TeV photon component compared to that in 
the sub-MeV component estimated for the TeV photon emitting GRBs claimed to 
have been detected so far. 

Before going into quantitative discussions in the following sections, it is necessary 
to specify the relativistic reference frames used in the calculations below. There 
are three reference frames of interest: (1) The “Wind Rest Frame” (WRF), which 
is comoving with the average ultrarelativistic “wind” shell, (2) the “lab” frame, in 
which the ultrarelativistic outflow moves with the average Lorentz factor T, and (3) 
the “observer’s frame” based on Earth. Quantities referred to in these frames will 
be indicated by the superscripts “WRF”, “lab” and “ob”, respectively. The “ob” 
frame quantities differ from the corresponding “lab” frame quantities only due to the 
cosmological redshift of the GRB source. Throughout this paper, physical quantities 
without any superscript for the reference frame will be understood to be referring to 
the WRF, unless explicitly stated otherwise. 


3 The Proton-Synchrotron Model of TeV Photon 
Production in GRBs 


3.1 The magnetic field strength in the wind rest frame 


We assume that there is a magnetic held ( B) within the GRB environment, the 
strength of which is determined by the condition that the total energy density in the 
system is approximately equipartitioned between the magnetic held energy density, 
B 2 /87 T, and total particle-radiation energy density, U. Following Totani [16, 17, 20, 32] 
and Ref. [12] we assume that U ~ U p ~ (m p /rn e ) LQ, where U 1 is the energy density of 
the low energy BATSE photon component (assumed to arise mainly through electron 
synchrotron radiation), and U p is the energy density in the form of protons. Note 
that all the above quantities refer to the WRF. We thus have 


B 2 

8tt 



( 1 ) 


where £b is the equipartition factor assumed to be of order unity. For a GRB at a 
given redshift z, the WRF energy density U 1 on the right hand side of Eq. (1) can be 
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expressed in terms of the observed variability time scale of the burst, Af ob , and the 
observed luminosity, E£ b , in the low energy (keV-MeV) BATSE component inferred 
from the measured fluence and the known distance (redshift), in the following way: 


rob _ 

L “ (1 + ^) 2 


r 4vr(r blb ) 2 r 2 d7 7 , 




( 2 ) 


where 

rjf b = 2T 2 cAf ob /(l + z) (3) 

is the lab-frame value of the characteristic “dissipation radius” where the internal 
shocks are formed and from where most of the radiation is emitted. Note that r'f b 
is the value of the dissipation radius as inferred by the observer at Earth from the 
measured variability time scale after correcting for the time dilation due to expansion 
of the Universe. 


Note also that E£ b as dehned above is the eAective 47r luminosity of the source assum¬ 
ing isotropic emission. If, instead, the source emits non-isotropically, for example, in 
a jet with opening solid angle U, then the actual luminosity would be a factor (Q/Att) 
smaller than E£ b . The same applies to the luminosity in the high energy component, 
Lh, used later. Obviously, our Anal results for the number of events per unit area 
in a detector on Earth and the associated signal-to-noise ratios are independent of 
whether or not the emission is isotropic, provided, of course, in the case of jetty 
emission, the emission jet from the source under consideration points toward Earth. 
Below we shall use the eAective 47r luminosity in all our calculations. 

Using equations (1) - (3), we can write 


B = 6.8 x 10 3 


d /2 i 

(rob N 

1^51, 

) V2 (1 + ^) 2 

rioo ( 

' A t° b j 
y sec ) 

1 


Gauss, 


( 4 ) 


where r 30 o = T/300, and T£ b 51 = /10 51 ergs 


-,-1 


3.2 Synchrotron photon spectrum 

The synchrotron radiation photons emitted by protons of energy E p in a magnetic 
held of strength B have the characteristic energy 


EL 


3 ehBE 2 p 
2 m^c 5 


(sin cc), 


( 5 ) 
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where e is the proton electric charge, and (sin a) = 7 t/4 is the average value of the sine 
of the pitch angle a of the proton. Using Eq. (4) for the magnetic held, we get the 
relation between the proton energy and its synchrotron photon energy in the WRF: 




( 6 ) 


where 


/ci = 1.3 x 10 8 


r 3/2 

1 300 _ 

#‘(£&i ) 1/4 U + *>' 


Note that we can use the general relation 


E ob = —E 
1 + z 


(7) 

( 8 ) 


between energy in the observer’s frame and and that in WRF, to express Eq. (6) in 
terms of observed energies: 


E° h ( E° h \ 1/2 

GeY _K 2 l v GevJ ’ 


(9) 


where 


k 2 = 2.3 x 10 £ 


r 2 

1 300 


' At ob \ 1 / 2 


sec / 


(b 4 {Li, 


51 


1/4 


(1 + z) 3 / 2 


( 10 ) 


/ r \i / 2 

Note that n 2 = (j—;) /ci- 

Now, let 


d ' n p a ( \ p 

dE p p VGeV / 


( 11 ) 


be the differential spectrum of the protons (number density of protons per unit en¬ 
ergy) accelerated in the internal shocks. The differential spectrum, dn 1 /dE 1 , of the 
synchrotron photons radiated by these protons is then related to the proton spectrum 
as 


fpJ,Ep)E p dE p — 


( 12 ) 


where f ps (E p ) is the fraction of energy of a proton of energy E p lost to synchrotron 
radiation, and E p and E 1 are related through Eq. (6). 
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The fraction f ps (E p ) is just the fractional energy loss of the proton of energy E p 
through synchrotron radiation during one expansion timescale of the wind, i.e., 


«£„) = 


t, 


exp 


where 


tp S (E p ) ’ 

1 (dE 


(E p 


tps(Ep) E p \ dt J pg 

is the inverse of the synchrotron loss time scale, with 

(^) P , (Bp) ~ C 9 (tv 5 ) 


(13) 


(14) 


E n 


B z 


(15) 


(assuming E p S> m p c 2 and isotropic distribution of the pitch angles), and 


1 r'fh 

f - 

texp - j, c 


(16) 


is the expansion time scale of the wind, both time scales being measured in the WRF. 

Since t ps (E p ) oc E “ 1 , we see that protons of energy above a “break” energy E p t, given 
by the condition 


tps(E p E pb ) t ex p, 


(17) 


will lose all their energy through synchrotron emission within one expansion time 
scale of the wind, thus giving 


fAEp) = | 

Using equations (14), (15) and (4) we get 


, if Ep < Epb , 

if E p > Ep b . 


(18) 


1 = 1.03 x l0-> ( L -) (1 +Z) E ^ 


which, together with equations (17), (16) and (3) gives 

( A t ob \ 

\ sec J 


\GeVj 


sec 


-l 


(19) 


Epb ~ 1.6 x 10 


r 5 

1 300 


(1 + z) 3 ^b (-^ L b 5 l ) 


GeV. 


( 20 ) 
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Equations (12), (11), (6), (7), (18) and (20) then give the synchrotron photon spec¬ 
trum in the WRF: 


where 


driy 

dE 7 



if E 1 < E ^ b , 
if E 1 > E lb , 



( 21 ) 


( 22 ) 

(23) 


is the “break” energy above which the synchrotron differential photon spectrum steep¬ 
ens by a power-law index 0.5 relative to the spectrum below that energy. 

The spectrum (21) cuts off above an energy E^ m&x given by 

JW = *r 2 (%=) 2 GeV, (24) 


where E p max is the maximum energy of the protons accelerated in the internal shocks. 
The acceleration time scale is f acc (E p ) rs-/ 27 TT]r L (E p )/c, where r L (E p ) = E p /(eB) is 
the Larmor radius and r) a factor of order unity. For reasonable ranges of values of 
the parameters of our interest, we have f acc < t exp < t ps for E p < E pb , f acc < t ps < 
fex|) for Epb < E p E Ep ^ rnax , and t ps ^ t acc for E p )> Ep r][ . !lx . (Here the equality 
signs in the time scale conditions refer to the corresponding equality signs in the 
energy conditions.) Thus, E pm . ax is in general larger than E pb , and is determined by 
synchrotron energy loss, i.e., by the condition 


^acc (E r ,max) — tps^Epy 


(25) 


Equation (25), together with equations (14), (15) and (4), gives 


B p> max 9.8 X 


8.n-V 2 - 


10 77 


r 3/2 / At ob \ 1/2 
1 300 


d /4 (i°L b 5 i) 


\ i/4 


( 1 +*) 


GeV, 


(26) 


and consequently, 

T 7 ,max = 54.73 t 7 _1 GeV. (27) 
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In the observer’s frame on Earth, equations (26), (27) and (8) give 

p5/2 / At °b\ V 2 

E° h max = 3 x 10 u r/- 1/2 - 300 ^ -GeV, (28) 

p d / 4 h&.) < 1 + 2 > 2 

and 

= 1-64 x 10 4 )f 1 (7S|) GeV . (29) 

In the numerical calculations below, we shall take rj — 1. Note that equation (28) 
shows that, depending on values of the various parameters involved, protons may in 
principle be accelerated to (observed) energies above 10 20 eV in GRBs. 


4 Optical Depth of High Energy Photons due to 
e + e~ pair Production and the Spectrum of Pho¬ 
tons at Source and at Earth 

The spectrum, eq. (21), of the high energy photons produced by the proton syn¬ 
chrotron mechanism, has to be corrected for optical depth of the produced photons 
within the source (“internal optical depth”) in order to calculate the spectrum of 
the high energy photons emerging from the GRB source. The emergent spectrum 
so calculated has to be then further corrected for optical depth of the emitted high 
energy photons in the intergalactic medium (“external optical depth”) on their way 
from the source to Earth, in order to calculate the spectrum of photons arriving at 
Earth. 


The dominant source of the optical depth of high energy photons within the GRB 
environment as well as in the intergalactic medium is e + e _ pair production due to 
collision of the high energy photons with other photons in the medium. A test photon 
of energy E t can produce an e + e” pair in a collision with a target photon of energy 
e greater than a certain threshold value eth given by 

2(m e c 2 ) 2 . 

e th = ~B~T\ -7n ) (30) 

E t ( 1 — cos 9) 

where 9 is the angle between the two photons’ directions of propagation. The cross 
section for this process for target photon energies e above the threshold is [33] 


u 77 (E t , e, 9) = —0t(1 - (3 ) 


(3 - /3 4 ) In 


1 + (3 


- 2/3(2 - (3 2 ) 


(31) 


where gt is the Thomson cross section and (3 = [1 — (eth/e)] 1 / 2 is the center-of-mass 
speed of the outgoing pair particles in units of c. 
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4.1 Internal Optical Depth and the Lab Frame Spectrum of 
Photons Emerging from the Source 

Let us first calculate the internal optical depth. The mean free path, £y 7 , of the test 
photon of energy E t within the GRB source can be written, in the WRF, assuming 
isotropic distribution of the photons in the WRF, as [34] 

1 r coo cIti 

C!(E?) = - / d(cos6>)(l - cos6>) / de — u 77 , (32) 

2 J J,cos 6) Cl(: 

where dn/de is the differential number density of the target photons (in the WRF) 
within the GRB source. 


From Eq. (30), we see that the lowest target photon threshold energy (which occurs 
for cos 6 = —1) in the WRF, for a given test photon of observed energy E° b = 
(T/(l + z))E t , is 


Gh ,min 


250 


/ 300 GeV\ 

l E? h ) 


^300 

1 + z 


eV, 


(33) 


which, in the observer’s frame takes the form 


e 


ob 

th,min 


/ 300 GeV\ 

l Ef ) 


r 3 oo \ 2 

1 + z) 


keV, 


(34) 


From the above two equations, it is clear that, for the high energy photons of our 
interest, the target photons for pair production within the GRB source are mainly 
those that have observed energies typically in the low (keV-MeV) energy BATSE 
range. We shall, therefore, assume that the dominant contribution to the internal 
optical depth due to pair production for the high energy photons of our interest 
comes from target photons that constitute the observed low energy BATSE photon 
component. With this assumption, then, we can write 


dn d<p oh 

de de oh 


(35) 


where d(f> oh /de oh is the differential spectrum of the observed BATSE photons, which, 
for most GRBs, is known to be well-fitted, within the BATSE energy range between 
e min ~ 20keV to e ° b aY ~ 2MeV, by a power-law with a break [35], 


d(P oh ^ f (ef) A 0h (e ob ) 31 , if e ob < ef, 
* ob * j (e ob ) _/3 \ if e ob > e% h , 


(36) 


where ~ 1 MeV is the break energy of the observed spectrum, and f3i ~ 1 and 

fdh ~ 2 . 
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The proportionality constant in (35) can be determined in terms of the observed 
luminosity L£ b defined in equation (2) through the condition 



f/ 7 


(1 + z? Lt ' 

47rc (r,if b ) T 2 


(37) 


together with the relation (8) between the energies in the WRF and observer’s frame. 


We can now perform the two integrals in equation (32) to calculate l~^. In doing 
this, it is convenient [34] to change variables from (e, cos#) to (e,s), with variable s 
defined as 


cE t (l — cos 6) 
2 (m e c 2 ) 2 


e 

— = s 0 z , 

^th 


(38) 


with 


s 0 = 


eEt 

(- m e c 2 ) 2 ’ 



(1 — cos 9). 


(39) 


Note that (3 = (1 — 1/s), so <v 77 now becomes a function of the single variable s, 
allowing us to write equation (32) as 


jr^ dee 2 ^[s°( e )] ’ ( 40 ) 

E t 

where 

rso(e) 1 (S /T 

</?[so(e)]= / sa(s)ds, with a(s) =-—. (41) 

J i ' 3 <tt 

The e and s integrals are now decoupled and are easily evaluated numerically. Finally, 
the internal optical depth T mt (E t ) for the high energy test photon of energy E t is given 
by the ratio of the wind expansion time scale t exp = r)j ab /cr and the mean two-photon 
collision time scale t 77 = ^ 77 /c, giving 


C(Et) = 


t 9 4 N 

m e c 


E t 


^lab 

Tint (Et) = -|r^ 77 • (42) 

The dependence of Ti nt on E t and on various other parameters in the problem can 
be easily extracted if we make the simplifying approximation <r 77 ~ (3/16 )ot = 
constant 3 . With this approximation, cx(s) = 1, and Lp [s 0 (e)] = (sq — l)/2, so the e 
integral in equation (40) can be performed analytically with the power-law form of 

3 Actually, for s 1, tr 77 oc (lns)/s, but the spectrum of target photons falls off rapidly for 
e > eth, so this approximation does not introduce large error. 
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dn/de given by equations (35), (36) and (37). For f3i — 1, f3h — 2, we get (in terms of 
quantities measured in the observer’s frame) 


rUEf) ~ 


where 


and 


3 (t t L^C 
647t c 2 At oh 


E? s = 


(1 + ^) 4 1 


Ef° 


T 6 3 m 2 c 4 ’ 


(1 + zf 

r 4 


c = 


4e? b 


+ 


12e? b 


'E%' 

E?\ 


+ 


2e' 


ob 


Z?ob 

lu\ 

&t,s 



l + z 


2 4 

m e c 

,-ob 


= 22.5 


+ In 


1 300 

TTz 


^.ob 


^ob 

e b 


1 MeV 

^ob 


GeV, 


for E° h < Ef° s 


for Ef° > £° b 


(43) 

(44) 

(45) 


is the energy below which, as equation (43) shows, r int oc E° h , while T int (E t ) “satu¬ 
rates” to a roughly constant value above this energy with two E° h dependent correc¬ 
tion terms one of which falls off as (E° h )~ 2 while the other increases with E° h only 
as In E° h . 


Equation (43) also shows the sensitive dependence of Ti nt on T. The results of our 
full numerical calculations of T; nt shown in Figures 1-4 clearly exhibit the expected 
dependence of Tint cm various parameters in the problem. 


With the internal optical depth calculated as above, the total number of high energy 
photons emitted by the GRB per unit time per unit energy as measured in the lab 
frame is given by 


dAd ab 

dE^dt™ 




(46) 


where dn 1 fdE 1 is as given by equation (21), and with E 1 = E]^°/Y. The total photon 
luminosity in the high energy component emitted by the GRB source (due to proton 
synchrotron radiation), as measured by a lab frame observer, is 


r lab — 
= 


dN lah 

E lab _ 7 _ dE 

7 dE^dt 1 ^ 7 


lab 


47T (7 


,lab\ 4 

d ) 


r 2 c 


E, 


dn~ 

dE, 


exp (-r int (E 7 )) dE. 


7 ? 


(47) 


where Ey tmax is given by equation (27). The exact value of E, hmin is not important 
as long as we take E ljmin < E lh since, the power-law index of the spectrum (21) 
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being less than 2 (for a p < 3 which we always assume to be the case), the dominant 
contribution to the energy integral in the region E j rnin < E y < E lb comes from 
its value at E lb . We use equation (47) to eliminate the constant A p appearing in 
equations (11) and (21) in terms of which we shall take as a free parameter. This 
completely fixes the lab frame emission spectrum of the high energy component given 
by equation (46). 


4.2 External Optical Depth and the Spectrum of High En¬ 
ergy Photons at Earth 

Due to expansion of the Universe, the energies of individual photons are decreased 
by a factor of (1 + z) and time intervals are stretched by a factor of (1 + z), where 
z is the redshift of the GRB source. Therefore, the spectrum of photons arriving at 
Earth (giving number of photons striking the top of the atmosphere per unit area per 
unit energy per unit time) is simply given by 


d,N° h 

dAdE° h dt oh 



1 dAd ab 

4t tD 2 , dEl^dt 1 ^ 

^ I 


(B“)exp (-r«,(£“)) , 


(48) 


where, on the right hand side, E l ^° = (1 + z)E° h , T ext (E l ^ h ) is the external opti¬ 
cal depth, clue to pair production through two-photon collision in the intergalactic 
medium, of a photon originating with energy E l ^ h at the source, and D z is the radial 
coordinate distance of the GRB source at redshift z. 


For a spatially flat Universe (which we shall assume to be the case) with rhv + f2 m = 1, 
where Da and Q m are, respectively, the contribution of the cosmological constant and 
matter to the energy density of the Universe in units of the critical energy density, 
3Hq/87iG, the radial coordinate distance D z is given by [36] 


= ( jr 

Mq 


dz' 


0 V^A + fhn(l + Z'Y 


(49) 


Here Hq is the Hubble constant in the present epoch. In our calculations, we shall 
use the currently popular values = 0.7, Q m = 0.3, and H 0 = 65 km sec” 1 Mpc” 1 . 

The main contribution to the external optical depth r ext of GeV-TeV photons comes 
from e + e” pair production due to collisions with the photons constituting the in¬ 
tergalactic infrared background, and has been studied by a number of authors. We 
shall use the optical depths given in [13], where a parametrization has been given for 
t(E 7 , z) for small z appropriate for our purpose. In our numerical calculations we 
use the optical depths given in Ref. [13] for the higher level of intergalactic infrared 


18 



background in order to have conservative estimate of the number of muons produced 
in the detector. 


5 TeV Photon Signal from GRBs in a ICECUBE 
Class Detector 

Photons of sufficiently high energy striking the earth’s atmosphere interact with the 
air nuclei to produce pions; these photo-pions are the major source of muon production 
by photons in atmospheric air-showers, among other possible physical processes like 
direct pair production of muons by photons or photo-prodnced charm decays. As first 
suggested by AH [24], it may be possible to detect the TeV photons from GRBs by 
detecting these (downward-going) muons created by the TeV photons in the Earth’s 
atmosphere in AMANDA/ICECUBE detectors. Onr aim in this paper is to explore 
this possibility within the context of a specific model of TeV gamma ray production 
in GRBs discussed in the previous sections. 


5.1 Muons in photon-induced air-showers and muon signal- 
to-noise ratio 


The number of muons with energy above E M in an air-shower produced by a single 
photon of energy E° h striking the top of the atmosphere can be parametrized [37, 24], 
for En in the range 100 GeV to 1 TeV, by the formula 


(50) 


where 6 is the zenith angle of the photon source, and all energies are in TeV units. 
The above parametrization is valid for E° b /E M >10 [24], 

The total number of muons with energy in excess of a threshold muon energy, E /l)th , 
in a detector of effective area A eff due to a single GRB of redshift z and observed 
duration T ob can then be written as 


rE° b 

N,(>E^ th ) = A eS T oh / dE ob - 

J E° h ., < 


d,N: 


ob 


E°\, 

7,th 


dAdE° h dt oh 


N,(E° b , >E,, th ) 


(51) 


where E° b h ~ 10 x E ^/cos 6 is the minimum photon energy needed to produce 
muons of energy E^the hi the atmosphere [24], 9 being the zenith angle of the photon 
source. 
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Equation (51) constitutes our signal, S = 7V /t (> E^ t h)- This must then be com¬ 
pared to the background of atmospheric muons created by cosmic rays in the Earth’s 
atmosphere. The differential spectrum of cosmic ray produced muons is given by [38] 


dN, 


/z,B 


0.14 E 


-2.7 


+ 


0.054 


dAdtdfldE u cm 2 ssr GeV 1 \ _i_ ( L1E ^ cos9 \ \ + ( l1 - e 'm cos6> 1 

1 f 1 ^ \ 115 GeV ) 1 ^ V 850 GeV ) 

The number of background muons in the detector is, therefore, 

dN^ B 


B = N^ b (>E^ th 


= A efl T oh AQ 


-'/Xjth 


dAdtdELdE, 


-dE , 


(52) 


(53) 


where we have included in the background the contribution only from the observed 
duration of the burst and from within a solid angle All = (7r<5d/180) 2 around the 
direction of the burst, 69 being the angular resolution (in degrees) of the detector. 
For AMANDA/ICECUBE type detector, 69 is typically a few degrees. We then define 
the signal-to-ratio, S/N, as 


S _ S-B 

N = 7b 


(54) 


5.2 Muon threshold energy and effective area of detector for 
downward-going muons 

For a muon to be detected by ICECUBE, it must have an energy > 100 GeV[21] 
within the detector volume. However, unlike in the case of the neutrino-induced 
muons which can be created by the neutrinos within or near the detector volume, the 
gamma ray induced muons in the atmosphere will have to have sufficient energy at the 
surface to reach the detector volume deep within ice (at a depth of > 1 km) and still 
have enough energy left to satisfy the trigger for the detector. This requires [39], for 
ICECUBE, a threshold muon energy of E ^ th > 250 GeV at the surface for a vertical 
muon. In this paper we take E ^ t h = 250 GeV in our numerical calculations of the 
muon signal-to-noise ratio for ICECUBE. 

We also need to specify what value to take for the effective area of the 
AMANDA/ICECUBE detectors for downward-going muons. As is well-known, these 
detectors are primarily designed to detect the upward-going muons created in the 
antarctic ice sheet by upward-going neutrinos from astrophysical sources, the main 
source of background for which is the upward-going muons due to upward-going at¬ 
mospheric neutrinos. Most of the PMTs constituting the optical modules (OMs) 
are, therefore, placed downward-facing, while only a small number of PMTs are kept 
upward-facing. The AMANDA-B10 detector [40], for example, has a total of 302 
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OMs, in 295 of which the PMTs are downward-facing, while 7 OMs have upward- 
facing PMTs. However, because of scattering of light in the ice surrounding the OMs, 
the OMs are not completely blind to light from its backward hemisphere. This effect 
makes each OM a relatively more isotropic light sensor in ice than in vacuum. For 
AMANDA-BIO detector, Ref. [40] estimates the effective relative sensitivity of an 
OM to be ~ 67% in the forward hemisphere and ~ 33% in the backward hemisphere. 
From this, we can estimate that the effective area of the AMANDA-BIO detector for 
downward-going muons would be ~ 1/2 of that for upward-going muons. 

The above discussion only serves to illustrate that although underice muon detectors 
like AMANDA/ICECUBE are designed to optimally detect the upward-going muons, 
the effective area of the detector for downward-going muons is still a significant frac¬ 
tion of its effective area for upward-going muons. 

Of course, the effective area for muon detection depends on the energy of the muon 
and also on the zenith angle ( 9 ) of the muon’s direction of propagation. Figure 23 of 
Ref. [40] shows that for the AMANDA-BIO detector, the effective area for upward- 
going (cosd ~ —1) muons is ~ 9000m 2 for muons in the energy range 100 - 1000 
GeV. The effective area then decreases with increasing cos 6 reaching a small value 
at cos 9 = 0 (horizontal direction). However, following the discussion in the previous 
paragraph, we expect that the effective area will increase again as we go above the 
horizon, i.e., as cosd increases from 0 toward 1 (directly downward-going), with the 
value of effective area roughly 1/2 of the corresponding value at the corresponding 
angle below the horizon. 

As will be clear from our results discussed in the next section, the effective area for 
AMANDA-B10, and for that matter even the AMANDA-11 detector, is not large 
enough to detect the TeV photon signal from GRBs with any reasonable degree of 
confidence within the context of the proton-synchrotron model discussed above, unless 
the GRB is extraordinarily bright in TeV photons. We shall, therefore, consider their 
detectability in the proposed ICECUBE detector which may reach an effective area 
of ~ 10 6 m 2 for upward-going muons (at cos 6 = —1). For the purpose of illustrating 
the kind of numbers for the signal-to-noise ratio we may expect in a ICECUBE class 
detector, we may assume that the ratio of the effective area for down-going muons 
to that for upward-going muons in ICECUBE will be roughly similar to that in the 
AMANDA-B10. We shall, therefore, take a value of ~ 5 x 10 5 m 2 as a fiducial value 
of the effective area for directly downward-going muons in ICECUBE for the muon 
energy range of our interest. Our results for the signal-to-noise ratio discussed in the 
next section can be easily scaled (as \JA eff ) for the actual effective area. 
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6 Results and Discussions 

We now discuss our main results. There are 11 free parameters in the model. These 
are: L^ 51 , £#, $, f3h, e£ b , Af ob , T ob , T, z, L\ f b , and a p . Out of these we fix L° b 51 = 1, 
£b = 0.5 , Pi — 1, Af ob = 0.5 sec , and T ob = (1 + 2 ) 10 sec , and study the dependence 
of our various results on the other remaining free parameters. Scalings of our results 
with respect to the fixed parameters are obvious from the relevant equations given in 
the earlier sections. 

Figures 1-4 show how the internal optical depth due to pair production within the 
GRB source as a function of the observed energy depends on various parameters that 
characterize the burst. The functional dependence on energy is as expected from 
the approximate analytic expression for T m given by equation (43), i.e., T m increases 
with energy up to a saturation energy given by equation (45) beyond which r in is 
approximately constant. 

The most sensitive dependence of r in is on T as indicated by equation (43) and as 
shown in Fig. 1: Larger values of Y give smaller values of r in , and thus depending on 
a moderately large value of Y in the range, say, 300-500, we can have r in < 1. 

Fig. 2 shows how Ti n (E° b ) depends on redshift of the burst. For a given value of 
the L£ b , higher redshift implies higher luminosity at the source. This in turn implies 
higher density of target photons for pair production inside the source, giving higher 
internal optical depth. Also, clearly, T m depends on the spectrum of the target photons 
within the source. We show the dependence of T m (E° h ) on two of the parameters that 
determine the target photon spectrum, namely, the break energy e£ b (Fig. 3) and the 
power-law index /3h of the spectrum beyond the break (Fig. 4). 

In Fig. 5 we display for illustration a typical predicted spectrum of the high energy 
component of a GRB due to proton-synchrotron process, where we also show the 
observed low energy “BATSE” spectrum for comparison. The effects of the internal 
and intergalactic optical depths due to pair production on the spectrum of the high 
energy component are also shown. 

In Figures 6-8 we show how the signal-to-noise ratio in the detector depends on 
various parameters that characterize the burst. In all these figures, as also in Fig. 9, we 
assume an effective area of the detector for downward-going muons to be ~ 5 x 10 5 rri 2 , 
assume the source to be at zero zenith angle, and take the threshold muon energy (at 
the surface) for detection by ICECUBE to be 250 GeV. Figure 6 shows the dependence 
of the signal-to-noise ratio on the redshift of the GRB for three different values of the 
emitted luminosity in the high energy component at the source, L}f b . Note that L^ b , 
as defined in equation (47), is the effective 47 t luminosity including the correction due 


22 



to internal optical depth. It is clear that for reasonable values of other parameters, 
detectability of the high energy component for GRBs at redshifts beyond about 0.05 
requires excessively high values of L \| b upward of ~ 10 58 ergs/ sec. Dependence of the 
signal-to-noise ratio on the Lorentz factor T of the underlying fireball model of the 
GRB is shown in Fig. 7. From Fig. 7, we see that TeV photons from a GRB at redshift 
z ~ 0.1 is detectable with a signal-to-noise ratio of ~ 10 if L]f b ~ 5 x 10 5 ' ergs/ sec 
and T ~ 400. 

In the proton-synchrotron model under consideration, the high energy photon spec¬ 
trum depends on the spectrum of the protons accelerated within the source. In Fig. 8, 
we show the dependence of the signal-to-noise ratio on the power-law index ( a p ) of 
the differential spectrum of the protons, for various values of the redshift of the GRB. 
As expected, a steeper proton spectrum (i.e., larger values of a p ) gives a steeper 
synchrotron photon spectrum, which gives smaller number of high energy photons 
emitted from the source and hence smaller signal-to-ratio in the detector. This de¬ 
pendence on a p is summarized in Fig. 9 where we show the behavior of the minimum 
luminosity in the high energy component, Lj| b , required for detection with a signal 
to noise ratio of 5 or larger, as a function of a p , for various values of the redshift z of 
the GRB. Clearly, as expected, a harder proton spectrum (i.e., smaller values of a p ) 
gives better prospect for detectability of the TeV photons in the proton-synchrotron 
model. 

From the above discussions while it is clear that, within the context of the proton 
synchrotron model, the energetic requirement for detectability of a possible high en¬ 
ergy ( xi TeV) photon component of GRBs in ICECUBE class detector is admittedly 
rather high, the required total energy estimates are not implausible, especially for 
reasonably close-by GRBs at redshifts < 0.05 or so, as seen from Fig. 9. 

We may compare the luminosity requirements discussed above for TeV photons 
with that for ultrahigh energy protons in the scenario where GRBs are sources of 
UHECR [18]. In order to explain the observed flux of UHECR, one needs a typi¬ 
cal GRB to emit a total energy of ~ few x 10 53 erg in UHE protons with energy 
> 10 19 eV, assuming that GRB rate evolves with redshift like the star formation rate 
in the Universe, with a local (z = 0) GRB rate of ~ 5 x 1CT 10 Mpc -3 yr —1 [41]. Note 
that the above number for the total energy in UHE protons refers to the total energy 
escaping from the GRB source in the form of UHE protons. The total energy in 
the UHE proton component produced within the GRB may be significantly higher, 
depending on the various energy loss processes of the UHE protons within the GRB 
source. If synchrotron radiation is the dominant energy loss process for the UHE 
protons within the GRB source, and if this process is very efficient, and further if the 
pair-production optical depth of the resulting GeV-TeV photons within the source is 
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sufficiently small, then the escaping GeV-TeV photon luminosity of the source may 
even be larger than the escaping UHE proton luminosity. 

An indicator of the efficiency of the proton-synchrotron process is the value of the 
energy E pbl given by equation (20), above which all protons lose all their energy 
through synchrotron radiation. Thus, smaller the value of this energy, higher is the 
synchrotron photon yield for a given spectrum of the protons. Specifically, if we 
assume that the internal pair-production optical depth of the high energy photons 
is negligible, then for a proton spectrum with index a p = 2, it is easy to show that 
the condition that the escaping high energy photon luminosity (due to synchrotron 
radiation of protons) be higher than the escaping proton luminosity is 

E pb <2.72 (£ p>min E p , max ) 1/2 , (55) 

where E pm - m and E Ptma x are the minimum and maximum energies of the protons, 
respectively. For appropriate choice of various parameters characterizing a GRB such 
as (measuring the internal magnetic held), T and so on, the above condition may 
well be satisfied. 

Lastly, we remind the reader that throughout the above discussion, the luminosity 
values we discussed refer to the isotropic effective 47 t luminosities; as already men¬ 
tioned, for collimated emission the required luminosity values would be appropriately 
smaller by a factor of Q/4n, where is the solid angle of the collimated jet. 

7 Summary and Conclusions 

It is possible that GRBs emit not only sub-MeV photons as detected in satellite- 
borne detectors, but also higher energy photons extending to TeV energies. For a 
GRB photon spectrum falling with energy, as is usually the case, the non-detection of 
TeV photons in satellite-borne detectors could be due to their limited size. However, 
as first suggested in Refs. [23, 24], TeV photons from GRBs might be detectable in the 
existing under-ice muon detectors such as AMANDA and future detectors such as the 
proposed ICECUBE, by detecting the muons in the atmospheric showers created by 
TeV photons. In this paper, we have made detailed calculations of the detectability 
of possible TeV photons from individual GRBs by AMANDA/ICECUBE detectors 
within the context of a specific model of such high energy gamma ray production 
within GRBs, namely, the proton-synchrotron model, which requires protons to be 
accelerated to ultrahigh energies > 10 20 eV within GRBs. In this model, the high 
energy component is distinct from, but may well be emitted in coincidence with, the 
usual “low” (keV-MeV) energy component observed by satellite-borne detectors. We 
have calculated the expected number of muons and the signal to noise ratio in these 
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detectors due to TeV gamma-rays from individual GRBs for various assumptions on 
their luminosity, distance (redshift), Lorentz Gamma factor of the underlying fireball 
model, and various spectral characteristics of the GRBs, including the effect of the 
absorption of TeV photons within the GRB as well as in the intergalactic infrared 
radiation background. The intergalactic absorption of TeV photons essentially pre¬ 
cludes detection of TeV photons in the currently operating AMANDA detector for 
any reasonable values of the luminosity in the high energy component, but they may 
well be detectable in the proposed ICECUBE detector which may have an effective 
area for downward-going muons a factor of 100 larger than that in AMANDA. How¬ 
ever, even in ICECUBE, only relatively close-by GRBs at redshifts < 0.05 or so can 
be expected to be detectable with any reasonable degree of confidence. The required 
isotropic luminosity of the high energy photon component is upward of 10 56 ergs/ sec 
and are generally found to be more than 3-4 orders of magnitude higher than typical 
estimated isotropic luminosities in the keV-MeV BATSE energy band. Such high 
luminosities may not, however, be impossible in some GRBs depending on various 
parameters that characterize the GRB. 
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Photon Energy in the Observer’s Frame (GeV) 


Figure 1: Internal optical depth as a function of the photon’s energy in the observer’s 
frame, for various values of T as indicated. Values taken for other relevant parameters 
are: L£ b 51 = 1, (3i — 1, (3h = 2.25, e£ b = 0.5 MeV, At oh = 0.5 sec, and z = 0.1. 
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Figure 2: Internal optical depth as a function of the photon’s energy in the observer’s 
frame, for various values of the redshift z of the GRB as indicated. Values taken for 
other relevant parameters are: L^ 5l — 1, /3 t = 1, / 3 h = 2.25, e£ b = 0.5 MeV, At ob = 
0.5 sec, and T = 300 . 
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Figure 3: Internal optical depth as a function of the photon’s energy in the observer’s 
frame, for various values of the break energy of the observed BATSE spectrum, e£ b , 
as indicated. Values taken for other relevant parameters are: L^ 51 = 1, $ = 1, (3h — 
2.25, Af ob = 0.5 sec, T = 300, and * = 0.1. 
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Photon Energy in the Observer’s Frame (GeV) 

Figure 4: Internal optical depth as a function of the photon’s energy in the observer’s 
frame, for various values of (3h as indicated. Values taken for other relevant parameters 
are: L^ 51 = 1, A = 1, eg b = 0.5 MeV, A t oh = 0.5 sec, T = 400, and z = 0.3. 
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Figure 5: A typical predicted differential spectrum of the high energy component 
of a GRB clue to proton-synchrotron process. The effects of the internal and inter- 
galactic optical depths due to 77 —> e + e _ process are indicated. The observed low 
energy “BATSE” spectrum is also shown for comparison. Values taken for the various 
relevant parameters are: L ° b 5{ — 1, £ B — 0.5, L \| b = 10 54 ergs/sec, a p = 2.2, = 

1, p h = 2.25, ef = 1 MeV, Af ob = 0.5 sec, T = 300, and z = 0.1. 
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Redshift of the GRB 

Figure 6: Signal to noise ratio in a ICECUBE class detector as a function of the 
redshift of the GRB for various values of the luminosity in the high energy component 
emitted from source, L^ b , as indicated, including the effects of the internal as well as 
external (intergalactic) optical depths. The zenith angle of the GRB source is assumed 
to be zero, and the effective area of the detector for downward-going muons is taken to 
be ~ 5 x 10 5 m 2 . The threshold muon energy (at the surface) is taken to be 250 GeV. 
Values taken for other relevant parameters are: L^i = 1, £b = 0.5, a p = 2.2, /3i = 
1, (3 h = 2.25, ef = 1 MeV, Af ob = 0.5 sec, T = 300, and T oh = (1 + z)10sec. 
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Lorentz factor 


Figure 7: Signal to noise ratio in a ICECUBE class detector as a function of the 
Lorentz factor T of the relativistic outflow characterizing the GRB, for various values 
of the luminosity in the high energy component emitted from source, L\ f b , as indi¬ 
cated, including the effects of the internal as well as external (intergalactic) optical 
depths. The zenith angle of the GRB source is assumed to be zero, and the effective 
area of the detector for downward-going muons is taken to be ~ 5 x 10 5 m 2 . The 
threshold muon energy (at the surface) is taken to be 250 GeV. Values taken for 
other relevant parameters are: LL b 51 = 1, = 0.5, a p = 2.2, f3i = 1, j3h = 2.25, e£ b = 

IMeV, A t oh = 0.5 sec, z = 0.1, and T oh = (1 + ^)10sec. 
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Differential Spectral Index of Protons 


Figure 8 : Signal to noise ratio in a ICECUBE class detector as a function of the 
power-law index (a p ) of the differential spectrum of protons accelerated within the 
GRB source, for various values of the redshift 2 of the GRB as indicated, including 
the effects of the internal as well as external (intergalactic) optical depths. The zenith 
angle of the GRB source is assumed to be zero, and the effective area of the detector 
for downward-going muons is taken to be ~ 5 x 10 5 m 2 . The threshold muon energy (at 
the surface) is taken to be 250 GeV. Values taken for other relevant parameters are: 


L£ b 51 = 1, A = 1, Ph = 2.25, ef = 1 MeV, = 0.5, A t oh = 0.5 sec, T = 300, Ljf = 
1 Q5S er g S / seC; an( ] = (1 -f ^)10 sec . 
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Figure 9: The minimum luminosity in the high energy component, L\ f b , required for 
detection with a signal to noise ratio of 5 or larger in a ICECUBE class detector, 
as a function of the power-law index ( a p ) of the differential spectrum of protons 
accelerated within the GRB source, for various values of the redshift z of the GRB 
as indicated, including the effects of the internal as well as external (intergalactic) 
optical depths. The zenith angle of the GRB source is assumed to be zero, and the 
effective area of the detector for downward-going muons is taken to be ~ 5 x 10 5 m 2 . 
The threshold muon energy (at the surface) is taken to be 250 GeV. Values taken 
for other relevant parameters are: L^ 51 — 1, (3 t — 1, (3h = 2.25, e£ b = 1 MeV, = 
0.5, A t oh = 0.5sec, T = 300, and T oh = (1 + z) 10sec. 
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